The electronic structure around As antisite near (110) surface of GaAs 
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The electronic structure around a single As antisite in GaAs is investigated in bulk and near 
the surface both in the stable and the metastable atomic configurations. The most characteristic 
electronic structures of As antisite is the existence of the localized p-orbitals extending from the As 
antisite. The major component of the highest occupied state on As antisite in the stable configura- 
tion is s-orbital connecting with neighboring As atoms with nodes whereas that in the metastable 
configuration is p-orbital connecting without nodes. Localized p-orbitals on the surrounding As 
atoms around the As antisite exist in every configuration of As antisite. Such features are retained 
except the case of the As antisite located just in the surface layer in which the midgap level is 
smeared into the conduction band and no localized states exist near the top of the valence band. 
Scanning tunneling microscopic images of defects observed in low-temperature grown GaAs, possi- 
bly assigned as As antisite, the origin of the metastability, and the peculiarity of the defects in the 
surface layer are discussed. 

PACS numbers: 73.20.Hb, 61.72.Ji, 68.37.Ef, 71.55.-i 
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I. INTRODUCTION 

Bulk GaAs crystals grown under As- rich conditions are 
known to contain As-related defects called EL2 centers 
1] that form a deep double donor level at the middle of 
the band gap. The EL2 centers compensate residual shal- 
low acceptors thus giving the crystals a semi-insulating 
property useful for device applications. The EL2 cen- 
ters induce a characteristic optical absorption extend- 
ing to the sub-band gap region with a spectral hump 
peaking around 1.2 eV. The most peculiar feature of the 
EL2 centers is the photoquenching effect in which the 
EL2-originated optical absorption becomes diminished 
when the crystal is optically illuminated with an intra- 
center excitation at low temperatures. The centers 
once quenched are, however, recovered if the crystal is 
heated or illuminated with a different infrared light. |3( 
This reversibility and the lost of the effect at high tem- 
peratures indicate that the EL2 centers, when they are 
electrically neutral, are transformable between the stable 
(EL2° ) and the metastable (EL2* ) state. 

For the atomic structure of EL2 centers, the most com- 
monly accepted model is isolated arsenic antisite defects 
(Asca ) ■ First principle calculations 0, |f| based on the 
Asca model showed that the metastable EL2* may be a 
close pair of an interstitial arsenic (Asj ) and a gallium 
vacancy (Vca ) that is formed by the displacement of the 
antisite As atom from the lattice point to a puckered 
interstitial position in a [111] direction leaving a Vca 
behind. The point symmetry lowering from Td to C3 V by 
the displacement was supported by optical absorption ex- 
periment. However, the atomic structure of EL2° and 
EL2* have been controversial since some defects, though 
exhibiting features of EL2 centers, were found to show no 
photoquenching effect. 0, The presence of such vari- 
ations in EL2-like centers arose an argument that the 
EL2 centers form a 'family' including various structures 
rather than have a unique atomic configuration. [t| 



Low-temperature grown (LT-) GaAs crystals homoepi- 
taxially grown under an excess As pressure contain a high 
density of As-related defects. Feenstra and his cowork- 
ers |lOj applied scanning tunneling microscopy (STM) 
to atomic level observations of individual point defects 
abundantly found in LT-GaAs epifilms. Four different de- 
fect contrasts observed in different sizes were attributed 
to Asca atoms located in different depths from the sam- 
ple surface. Recently, some of the present authors [llj 
provided direct evidence for the defect being EL2: They 
found that the STM contrasts of the defects in LT-GaAs 
samples kept at a low temperature (90K) are drastically 
changed by light illumination with an excitation spec- 
trum nearly identical to that for the photoquenching ef- 
fect of EL2 centers in bulk crystal. The local density 
of states (£DOS) measured by scanning tunneling spec- 
troscopy (STS) at the defect sites shows that the donor 
gap state present in the normal state disappears in the 
metastable state, conforming the characteristics of EL2. 

Theoretically simulated STM images of Asc a defects 
show an apparently good agreement with experimental 
images as far as the defects in the second largest con- 
trast ^3 are concerned. However, no systematic 
comparison has been made for the STM image contrasts 
in other sizes and those in the metastable state, ^ij 

In this paper, we will report results of first principle 
electronic structure calculations of an Asca in GaAs bulk 
crystal and those near the surface, both in the stable and 
the metastable state. The main aim of the present study 
is to answer a naive question if defects near the surface 
that can be probed by STM may be substantially differ- 
ent or not in the structures and physical properties from 
the defects in the bulk crystals. In Sec. |Hj we describe 
the methods of calculations based on the atomic structure 
models. In Sec. IIIII we show the electronic structures of 
the highest occupied level in real space and discuss the 
nature of the level which differs depending on the depth 
of the Asca from the surface. Section IIVI is devoted to 
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a comparison between the simulated STM images and 
experimental ones and discussions about the effects of 
surface on the properties of Asc a centers in light of the 
calculated electronic structures. In Sec. we conclude 
this paper. 



II. COMPUTATIONAL DETAILS 

Since a considerable lattice expansion is present around 
Asca defects due to the antibonding nature of the dou- 
ble donor level in the band gap, |l5j it is needed to 
calculate the electronic structure allowing lattice relax- 
ation with sufficient accuracy. For this purpose, we 
used the density-functional (DFT) theory with the local- 
density approximation (LDA) employing two different 
methods for electronic structure calculations. One is 
the method using the norm-conserving pseudopotential 
of the Troullicr-Martin type with the plane wave basis 
set ^(J and the other is the tight-binding Linear Muffin- 
Tin Orbital (LMTO) method. H3 

In the case of the norm-conserving pseudopotential of 
the Troullier-Martin type, we adopt a supercell and the 
cut-off energy for the plane wave basis was set to be 
10 Rydberg. The k-point sampling in the Brillouin zone 
was done only at the T point. We relaxed atomic posi- 
tions, and calculated the electron density with the plane 
wave basis. 

The LMTO basis is localized around each atom and 
hence we can analyze the contribution of the atomic or- 
bital to a specific energy state in concern. The LMTO 
basis set is minimal and each basis function can be de- 
scribed by a small number of parameters so that we can 
perform calculations with relatively small memory and 
computing time. 

The procedure for calculating the electronic structures 
in the model system that is bounded with a surface and 
contains an AsQa defect is as follows: 

(1) Relaxing a relatively small lattice system with an 
Asca and a surface by the pseudopotential method. 

(2) Embedding the relaxed small system in a larger per- 
fect lattice with the atomic configuration around the 
Asca unchanged. 

(3) Calculating the electron density with T point in a su- 
percell by the pseudopotential method, and other proper- 
ties such as the E — k relation, the total density of states 
(DOS) and local density of states (£DOS) of each atom 
by the tight-binding LMTO method. 

For Asca defects in bulk crystals, we introduced an 
AsQa defect by replacing a Ga atom with an As atom in 
a 2 x 2 x 2 cubic unit cells containing 64 atoms. These 
atoms were then subject to relaxation by means of the 
molecular dynamics with plane wave bases. The lattice 
constant was fixed to the experimental value 5.654 A dur- 
ing the relaxation. We found that the distances between 
the AsQa and the neighboring atoms after relaxation are 
similar to the previous results. [THl Il8l| 

For calculations of the defect structure in the 
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FIG. 1: The definition of (a) the layer number and (b) three 
inequivalent directions to which the AsGa is displaced to four 
possible metastable positions. 



metastable state, we first displaced the Asc a atom ten- 
tatively by 1.2 A from the lattice point to a puckered 
interstitial position along a bond breaking [111] direc- 
tion, and then relaxed the lattice by means of the molec- 
ular dynamics with plane wave bases. The resultant dis- 
tance after the relaxation between Asi and a neighboring 
As atom on the axis of the symmetry Cjvwas 1.37 A, 
a little larger than the previous values. 0, 0, 0] The 
electron density were calculated for the perfect lattice of 
this size (64 atoms) and for lattices in an extended size 
of 4 x 2v2 x 3/\/2 (192 atoms) in which we embedded 
a relaxed 2x2x2 supercell containing an AsGa or an 
As; -Vca pair. In these calculations, we fixed the atoms 
beyond the third neighbors at the perfect lattice posi- 
tions. 

For Asca defects near a (110) surface, we prepared a 
slab lattice parallel to the (110) plane in which one sur- 
face was relaxed and the other surface was terminated 
with hydrogen atoms. For convenience, we refer to the 
surface layer without hydrogen, though buckled as shown 
in Fig.^a), as 'layer 1', the one layer down as 'layer 2', 
one more layer down as 'layer 3' and so on. The electronic 
structures of the system bounded with surfaces were cal- 
culated for various configurations: without Asca defects, 
with an Asca defect (the stable state) located in the layer 
1 to the layer 4. 

The procedure of lattice relaxation and electronic 
structure calculations are as follows. We first relaxed 
the surface bounded slab of a lateral size of 3 x 3/V2 
and a thickness of 7 layers including the layer of hydro- 
gen atoms. The lattice constant was fixed in this case 
too. The positions of atoms in the relaxed configuration 
displaced from those of the ideal atomic positions are 
smaller than 0.2 A, 0.1 A and 0.05 A for the first, the 
second and the third neighbor atoms of Asca , respec- 
tively. Some of the neighbor atoms were displaced con- 
siderably, especially near the surface. For example, for an 
Asca located on the layer 2, the third neighbor atoms on 
the layer 1 were displaced by 0.063 A, whereas the third 
neighbors on the layer 3 only by 0.019 A. Therefore, the 
atoms at shallow depths from the surface were not fixed 
tightly while the atoms further than third neighbors were 
fixed since all these atoms, even if relaxation was allowed, 
were found to be displaced by only a magnitude smaller 
than 0.07 A. After relaxing the small system, we em- 
bedded it in a larger system with a surface of a size of 
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4 x 4/V2 and 5 layers including the hydrogen layer. For 
an Asca on the layer 4, we embedded the relaxed system 
in a still larger system of the same surface size 4 x 4/V2 
but of 6 layers in thickness. 

For calculations of metastable Asi-Vca pairs, we 
should note that there are three inequivalent configura- 
tions regarding the displacement direction relative to the 
surface as shown in Fig.QJb). We hereafter call the "up" 
configuration for the displacement direction ascending to 
the surface, the "side" configuration for the direction par- 
allel to the surface, and the "down" configuration for the 
direction descending from the surface. The lattice re- 
laxed to a metastable configuration similar to that in the 
bulk crystal except for Asc a atoms on the layer 1, where 
the similar displacement of the surface Asc a atoms would 
destabilize the surface buckling and therefore is energet- 
ically unfavorable. 

According to Tersoff-Hamann theory, [2(| filled-state 
STM images under a negative sample bias of —2.3 V, 
the bias condition experimentally employed in the pre- 
vious study, were simulated by taking into account 
only the highest occupied state that should give the main 
contribution to the contrast. The image was calculated 
as a two-dimensional contour map of the iso-surface of 
electron density of the highest occupied state in vacuum 
out of the surface. 



III. ELECTRONIC STRUCTURE 
A. As antisite in bulk crystal 

1. Stable configuration 

Figures |3[a) and |2ft>) show the E — k relations calcu- 
lated by the LMTO method in a bulk system with a unit 
cell of 64 atoms. Figure|2fa) is that of the perfect lattice 
and 01b) that of the lattice containing an Asc a atom in 
the stable configuration. The band structure of the per- 
fect lattice is characterized by three bands with energy 
gaps, the As s band below —10.0 eV, the valence band 
between 0.0 eV and —7.5 eV, and the conduction band 
above the energy gap. An AsQa atom in the stable con- 
figuration introduces an occupied band with relatively 
small dispersion in the band gap as shown in Fig. 0{b). 
The DOS of perfect lattice and that of the lattice con- 
taining an AsQa in the stable configuration are shown in 
Figs. |2Id) and [2Je) respectively. Comparison between 
the two figures indicates that Asc a atom introduces the 
gap state mainly localized at the Asca atom and neigh- 
boring four As atoms. 

Figure [21 shows the iso-surface of electron density of 
the highest occupied state in the gap calculated by the 
plane wave basis for a system of 192 atoms containing an 
Asca in. the stable configuration. One may see the Td 
symmetric structure extending first along four [111] di- 
rections and then each further branching out along three 
[110] directions. The decreasing electron density with 



distance from the Asca shows the localized nature of the 
corresponding £DOS at the Asc a atom (Fig. Etc)). The 
iso-surface of electron density is shown also in Figs. 0Ja) 
and Bib). We note that the electron density at the 
Asca site is s-like (colored in red) and those on the near- 
est neighbor As sites are p-like (colored in blue) with the 
lobes heading toward the Asc a atom. These p-like iso- 
surface are respectively linked to the central s-like iso- 
surface with a distinct node. The electron density on the 
third nearest neighbor As atoms are also p-like (colored 
in green) with the lobes aligning in the same directions as 
those of the first neighbor As atoms and the fifth neigh- 
bors (colored in yellow) again forming nodes between 
them. Thus, the electronic structure of the highest occu- 
pied level in the stable configuration has a radial pattern 
that are characterized by the mutually anti-bonding As 
p-orbitals. 

The defect formation energy SI of Asc a in the unit cell 
of 64 atoms, is written in As-rich limit as |2lJ 

fl = E B - N c fi c - A^GaMGaAs ~ (^As ~ ^GaWsOulk) , (1) 

where Er> is the total energy of the system with AsQa , Me 
is the chemical potential of electron, /iGaAs is the energy 
per atomic pair of bulk GaAs, and /^As(buik) is the energy 
per atom of pure bulk As. N e (= 0) is the charge of the 
defect, and Aca and N\ s are the numbers of atoms of 
each species in the system. The defect formation energy 
is estimated to be 0.6 eV by using the T point sampling 
for the total energy calculation. However, this total en- 
ergy is underestimated since, as seen in the Fig. |2fb), 
r point is the bottom of the dispersion of the midgap 
impurity level. Since the peak of the midgap impurity 
level locates at 0.31 eV above its bottom, the amount of 
underestimated energy would be 0.6 eV. By adding it to 
the total energy, we estimate the defect formation energy 
to be 1.2 eV. This value is consistent with those of the 
previous calculations 1.4-1.8 eV. [101112 

2. Metastable configuration 

Figures |2fc) and|2If) show the E — k relation, the to- 
tal DOS and the ffiOS calculated by the LMTO method 
for a system with a unit cell of 64 atoms containing an 
AsQa atom in the metastable configuration (Asj -Vc a ) 
that was fully relaxed by the pseudopotential method. 
The analysis of the £DOS shows that the basis function 
localized at the Asi atom constitutes the main contribu- 
tion to the lowest unoccupied states introduced in the 
band gap between the valence band and the conduction 
band. Instead, there are no occupied states associated 
with localized orbitals at the Asi atom near the Fermi 
level. This implies that the lowest unoccupied states and 
the highest occupied states in the stable configuration 
exchange their energetic order in the metastable configu- 
ration. However, the position of the unoccupied level in 
the gap was found to be sensitive to the computational 
detail. For instance, a pseudopotential calculation with 
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FIG. 2: The band structure, the total density of states (DOS), and the local density of states (£DOS) in the energy range of 
the band gap in a relaxed system with a unit cell of 64 atoms. The energy zero is set to be at the position of the top of the 
bulk valence bands, (a) The E — k relation without defect, (b) The E — k relation with AsGa in the stable configuration, (c) 
The E — k relation with Asi-VGapair in the metastable configuration, (d) The DOSs without defect corresponding to (a). The 
solid line: total DOS, broken line: ffiOS of Ga, dotted line: ffiOS of As. (e) The DOSs with AsGa corresponding to (b). The 
solid line: total DOS, broken line: £DOS of AsGa, dotted line: averaged £DOS of neighboring four As atoms, (f) The DOSs 
with Asi-VGapair corresponding to (c). The solid line: total DOS, broken line: £DOS of Asi, dotted line: averaged £DOS of 
neighboring four As atoms. In the E — k relations, the dashed horizontal lines represent the Fermi energies. In comparison 
with the stable configuration (e), the Fermi energy Ef shifts in the metastable configuration (f) and the gap state becomes 
unoccupied. 



a larger supercell shifts this unoccupied level to a slightly 
higher energy position in the conduction band region. 

FiguresEIa) and[S{b) show the electron density of the 
corresponding highest occupied state of the system of 192 
atoms containing a relaxed Asj -Vca pair. The electron 
density of this state extends over a wider range and dif- 
fers substantially from the symmetric radial pattern in 
the stable configuration shown in Figs. E{a) and^b). 
The less localized feature is consistent with the large dis- 
persion of the highest occupied band that merges into 
the bulk valence bands as seen in Fig. [2fc). We should 
note that the iso-surface of the electron density at the 
Asi atom becomes p-like in the metastable configuration 
in contrast to the s-like feature in the stable configura- 
tion. The p-like iso-surface at the Asi is in parallel to 
those of the three neighboring As atoms with no nodes 
between them. The radial pattern characteristics of the 
stable configuration is broken to C3 V • Among the iso- 
surface branches extending in four [111] directions, the 
branch remained is only the one in the direction opposite 



to the displacement direction of the Asc a atom. The lobe 
of the branch, however, preserves the p-like anti-bonding 
character as in the stable configuration. 



B. As antisite located near (110) surface 

1. Stable configuration 

For slab systems with a surface, the electronic en- 
ergy levels were calculated only by the pseudopoten- 
tial method with plane wave basis. Except for the 
Asca located in the surface layer (the layer 1), the high- 
est occupied state remains amid in the band gap as in 
the bulk crystal. The absence of the gap state associated 
with Asca located in the layer 1 is shown also by the 
£DOS calculated by the LMTO for a system of 60 atoms 
with the surface size of 2 x 3/V2 and the thickness of 5 
layers terminated with hydrogens on one side. Figures 
0Jc)-(j) show the iso-surfaces of electron density of the 
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highest occupied states in the gap for an AsQ a located in 
the stable configuration in the layer 1-4. In all the cases, 
the features of the iso-surface is essentially the same as 
those in the bulk crystal (Figs. Efa) and^Jb)): the s-like 
feature at the Asc a site and the p-like ones at the neigh- 
boring As atoms. The radial pattern characteristics of 
the highest occupied state in the bulk crystal look as if 
it were just terminated with the surface. Even in the 
Asca located in the layer 2 (Figs. QJg) andQJh)), though 
the electron density is concentrated mainly on the first 
and the third neighbor As atoms near the surface with 
more distinct p-like features than in the bulk lattice, the 
atomic configuration around the Asca almost remains the 
same as in the bulk lattice. The Asca located in the layer 
3 (Figs. E{e) andQJf)), however, shows a subtle differ- 
ence in the iso-surface from that in the bulk crystal. The 
p-like iso-surface at the first nearest As atom extends 
considerably to Ga atoms in the layer 1 forming 7r-like 
bonds without a node. This may be an effect of the sur- 
face buckling: For the surface Ga, As in the layer 2, and 
Asca , the angle of Ga-As-Ascais at 87°, which is ap- 
proximately right angle, so that the p-orbitals on the As 
atoms can easily form tt bonds. The situation is quite 
different when the As antisite is located in the layer 1, 
for which Figs. 0Ji) and0]j) show the iso-surface of the 
electron density. The cause of such a surface effect on 
the electronic structure is discussed in Sec. IIVI 



"down" configuration gives a better agreement of the cal- 
culated STM image with that of experiments. Therefore, 
only the "down" configuratin is discussed for a possi- 
ble metastable state of As; -Vc a pairs in this section. 
When the As antisite is located just on the surface, no 
metastable configuration is found in simulation. 

Figures Efc)-(f) show the iso-surfaces of the electron 
density of the highest occupied state when the Asi -Vca 
pair is formed in the layer 3 and the layer 2. It is notable 
that the electron density increases exclusively along the 
direction opposite to the displacement direction of the 
Asca ; which is essentially the same as in the bulk lattice. 
The iso-surface localized on the Asj and the neighboring 
three As atoms in the puckered configuration has com- 
monly p-like lobes bonding with each other without nodes 
as in the bulk lattice. However, the degree of localiza- 
tion of the highest occupied state is slightly dependent on 
the depth of the Asi atom. When the Asi atom is located 
in the layer 2 (Figs.OJe) and|S{f)), the highest occupied 
state is well localized: The electron density is distributed 
mainly along the [llT] direction. When the Asi is located 
in the layer 3 (Figs.UJc) and EJd)), the electron density 
is less localized having a distribution over atoms outside 
the branches. Except for these fine details, the crystal 
surface has little effect on the electronic structure of the 
metastable configuration as in the stable one as far as 
the Asi atom is located below the layer 1. 



2. Metastable configuration 

As mentioned in Sec. [HI there are three inequivalent di- 
rections to which the antisite As atom could be displaced 
to form a Asi -Vc a pair. Since the STM experiment [Tl| 
shows that the defect contrasts in the metastable state 
are all symmetric with respect to the [001] axis, the "up" 
or "down" configurations could be a candidate for the 
metastable state. As shown later in Sec. IIVI only the 




FIG. 3: The typical iso-surface of the electron density of the 
highest occupied band in the stable configuration in bulk. 



IV. DISCUSSION 

A. VALIDITY OF EMBEDDED SYSTEMS 

In constructing the atomic configuration, two approx- 
imations are employed. First, we relaxed only the atoms 
in smaller supercell and fixed the lattice constant. The 
expansion of the lattice constant by the introduce of 
Asca is 1 %. 23] It is relatively small in comparison 
with the change of the distance 6 % between Asca and 
its nearest neighbor. The atomic configuration around As 
interstitial is well converged in the system of 65 atoms, 
[23| and, therefore, we presume that it is also the case 
in Asca ■ Although we actually tried to calculate the 
electronic structure without the relaxation of the atoms, 
there remains the features of the electronic structures 
discussed in Figs. 4 and 5. From these considerations, 
we conclude that little error is introduced by fixing the 
lattice constant. 

One may consider that the calculated wavefunction 
could be sensitive to the number of layers. We calculated 
the electronic structure of the system with smaller sur- 
face area and more layers, and found that the essential 
characteristics of wavefunctions are already converged. 
Therefore, the electronic structure is not sensitive to the 
number of layer, or the surface size. 



6 



(a) ioa 




FIG. 4: The iso-surface of the electron density of the highest occupied level of AsGa in the stable configuration calculated by 
pseudopotential method. The figures in left and right column are the top and side view of the electronic density, respectively. 
Figures (a) and (b) represent the iso-surfaces of electron density 1.37 x 10 -3 e/A 3 in the bulk crystal. Black spheres represent 
As atoms, and white spheres Ga atoms. The color red, blue, green and yellow, represent the wavefunction on the central 
AsGa , the first, the third and the fifth neighbor As atoms, respectively. The wavefunction has little amplitude on the Ga site. 
The average of the electron density of this level is 2.31 x 10~ 4 e/A 3 . In the stable configuration, there are nodes between 
the s-orbital of the AsGa (red) and p-orbitals of neighboring four As atoms. Figures (c) to (j) show the iso-surface of electron 
density 2.26 x 10~ 3 e/A 3 for an As Ga in the layer 4 ((c), (d)), the layer 3 ((e), (f)), the layer 2 ((g), (h)), and the layer 1 ((i), 
(j)). The average electron density of this level is 1.73 x 10 -4 e/A 3 . 
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FIG. 5: The iso-surface of the electron density when Asi is located at the metastable position in the bulk crystal ((a), (b)), 
in the third ((c), (d)) and the second layer ((e), (f)). The figures in left and right column are the top and side view of the 
electronic density, respectively, as Fig. 2] In the metastable configuration in bulk crystal ((a), (b)), there are no nodes between 
the p-orbital of Asi (colored in red) and p-orbitals (colored in gray) of neighboring three As atoms. See also the caption of Fig. 
4. 



B. SIMULATED STM IMAGES 

1. STM image in the stable configuration 

Figures EJa) , HJb) , and EJc) show the simulated filled- 
state STM images of stable Asc a atoms in different 
depths from the surface. The images were calculated by 
assuming a negative sample bias and hence considering 
only the iso-surface of the highest occupied state. The 
image of the AsGa m the layer 2 is in good agreement 
with the previous results. [TJ As mentioned in Sec. II I II 
even when the Asc a is located near the surface, the wave- 
function (iso-surface) of the highest occupied state does 
not differ much from that in the bulk lattice with the 
characteristic radial pattern. Owing to this fact, when 
the As Ga is located in the layer 2 (Fig. Eta)), the STM 
contrast of the defect core arises from the wavefunction 
concentrated at As atoms in the surface layer (the layer 
1) that are neighboring to the Asca > and the satellite 
contrasts near the core from the wavefunction at the 
third As neighbors also in the surface layer (the layer 
1). When the AsQ a is located in the layer 3 (Fig. EJb)), 
the highest occupied state has a large density on the sur- 
face just above the Asc a atom and the fifth As neighbors 



in the surface layer, giving the satellite features. For the 
Asca located in the layer 4 (Fig. 0c)), the satellite peaks 
arise from the seventh As neighbors. Thus, since the 
satellite peaks originate in the radial pattern of the iso- 
surface of the highest occupied state, they become more 
remote from the core contrast with increasing depth of 
the Asca from the surface. 

Feenstra et al. argued (later referred to as FWP) 
that the defect in the STM contrast in their largest size 
(type A) should be assigned to an Ascain the layer 1, 
the defect in the second size (type B) to one in the layer 
2, the third size (type C) to one in the layer 3, and the 
fourth size (type D) to one in the layer 4. Their criterion 
used for judging the depth of the defect was the relative 
position of surface As atom contrasts with respect to the 
defect core contrast. Geometrically, if the Asc a atom is 
located in an odd-numbered layer, the mirror symmet- 
ric (110) plane passing the AsGaatom must cut halfway 
between As atoms in the surface layer. In contrast, if 
the Ascais located in an even-numbered layer, the mir- 
ror plane must pass an As atom in the surface layer. The 
assignment by FWP of type B defect as an Asg a in the 
layer 2 agrees with the assignment by Capaz et al. [IJ As 
for an Ascain the layer 1, Ebert et al. reported a defect 



image that is quite different from type A defects by FWP, 
though such defects were not observed in our experiments 
probably due to the excessively large magnitude of the 
bias voltage. [24| Figures HJd) andHJe) show experimen- 
tal STM images typically observed at a sample bias of 
—2.3 V. Our experimental images shown in Figs. EJd) 
andH^e) seem identical to respectively the type B and C 
defects in experiments by FWP. [l0| 

Another criterion for depth assignment is the distance 
between the two satellites from the core contrast as stated 
above. The order of the experimental images in Fig. 
are tentatively arranged according to this criterion. Al- 
though the quantitative agreement is not perfect, we 
could assign the image shown in Fig. |5fe) to an Asc a in 
the layer 3, and image Fig.^Jd) to an AsGa in the second 
layer. The assignment of the type B defect (Fig. Efd)), 
the same as proposed by Capaz et al., |3| seems rea- 
sonable because the agreement of the defect contrast in 
the metastable state is better than otherwise as shown in 
what follows. 



2. STM image in the metastable configuration 

Figure [7Jc) shows the experimental STM image of the 
defect in the metastable state that underwent a change 
from the contrast shown in Fig. EJe). The image in the 
metastable state is characterized by the diminish of the 
satellite peaks and the concomitant appearance of a new 
contrast elongated along the surface As atom along [110] 
row on the side opposite to the diminished satellites. As 
shown in Sec. IIIII the branches of the radial wavefunction 
of the highest occupied state disappear upon transforma- 
tion to the metastable configuration, except for only a 
branch extending in the [111] direction opposite to the 
displacement of the Asc a atoms. This means that if the 
Asca atom near the surface was displaced to a "side" po- 
sition, the STM image would lose the mirror symmetry 
about the (110) plane. If the Asc a atom was displaced to 
the "up" position, the peak position move in the oppo- 
site direction to the one experimentally observed. Both 
of these contradict the experiments. [llj 

This conclusion is different from the calculations by 
Zhang, 0] where the "side" configuration in the sec- 
ond layer is the most stable energetically among three 
possible configurations of the metastable state. Our cal- 
culation shows the similar energy barrier in the path of 
Asca from the metastable to the stable positions in these 
configurations. The energy barrier heights are 0.5 eV in 
the "side" configuration in the second layer, 0.3 eV in 
the "down" configuration in the second layer (Fig. da)), 
and 0.3 eV in the "down" configuration in the third layer 
(Fig-Elk))- Therefore, it is still a open question why only 
one metastable defect image has been observed in STM 
experiment that is not energetically favorable. 

Figures[7Ja) and[7fb) show the simulated STM images 
of As; -Vca pairs located in the layer 2 and in the layer 
3, respectively, with the As; in the "down" configura- 
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FIG. 6: The iso-surface of the electron density 8.32 x 
10 _4 e/A 3 of the highest occupied band when AsGa locates 
near the crystal surface. The AsGa locates in (a) the layer 2, 
(b) the layer 3, and (c) the layer 4. The color becomes darker 
with the iso-surface approaching the vacuum. It can be seen 
that the distance between two satellite peaks becomes larger 
as AsGa locates at deeper positions. Figures (d) and (e) are 
experimental STM images. 



(a) 




FIG. 7: The iso-surface of the electron density of the highest 
occupied band when Asi locates at metastable positions in 
(a) the layer 2 and (b) the layer 3. The satellite peaks in the 
stable configuration vanish (Compare with Fig. EJ. Figure 
(c) is an experimentally observed STM image, of defects as 
Fig. Ele) in the metastable state. 



tion. The agreement between simulation and experiment 
is better if the Asc a is assumed to be located in the layer 
3 as tentatively assigned above. 
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C. Origin of metastability 

The physical origin of the metastability of Asi-Vca 
pairs can be understood if one examines the change in the 
bonding character upon the stable to metastable trans- 
formation. As shown in Sec. IIIII the Asca configuration, 
though lower in total energy, bears an inherent instabil- 
ity due to the anti-bonding nature of the chemical bonds 
with the surrounding atoms. The global expansion of the 
lattice by the presence of Asc a defects is, thus, owing to 
this anti-bonding nature of the electronic structure as- 
sociated with the defects. [23J On transformation to the 
metastable state, the electronic energy of the s-orbital 
associated with the Asc a atom that is anti-bonded with 
the nearest As atoms is lifted and becomes emptied. In 
its place, a gap state associated with an atomic orbital of 
p-character originating in the displaced As atom is drawn 
from the conduction band, and it becomes occupied to 
form bonds with three neighboring As atoms. Thus, some 
of the increase of the total energy on transformation is 
canceled by the newly formed chemical bonds which con- 
tributes to the affinity of the interstitial As atom with 
the surrounding. 

Since the energy gain by such bond reconstruction 
should become larger with decreasing distance between 
the atoms, the change of the dominant orbital of the 
Asca from s to p-orbital, and the disappearance of the 
anti-bonding node between Asi and the surrounding As 
atoms approached by Asca are considered to be the rea- 
sons why the metastable state exists. 

D. As antisite in (110) surface 

As shown in Figs. 01 and the presence of surface has 
no significant effects on the electronic structures associ- 
ated with AsGa defects as far as they are located deeper 
than the layer 2. It is only when the Asc a is located in 
the surface layer (the layer 1) that no level is formed in 
the band gap. Due to the absence of the gap state or the 
spatial spread of the occupied states, AsQ a defects in the 
layer 1 give rise to no localized STM contrast at negative 
sample biases, in agreement with the results obtained by 
Ebert et al. 

The surface buckling [25| which is reproduced in the 
present calculation as well is induced by the electron 
transfer in the Ga dangling bond which is higher in en- 
ergy to the As dangling bond. This causes a change of the 
orbital character from sp 3 in the perfect lattice to sp 2 -like 
on the surface Ga atom and to p 3 -like on the surface As 
atom, thereby resulting in the retraction of the Ga atom 
and the protrusion of the As atom from the surface. Con- 
comitantly, the As dangling bond state smears into the 
valence band, and the Ga dangling bond state into the 
conduction bands. It is expected that when Asca defect 
is located in the surface (layer 1), Asca reverts to the un- 
buckling position as noticed in a figure of the paper by 
Ebert et al. 



To test this idea, we calculated a model system of 60 
atoms with the surface size 2 x 3/y/2 and with 5 layers 
including the hydrogen layer by the LMTO method. This 
model had a clean (110) surface on one end, and the 
other surface was terminated by hydrogens. Then, AsQa 
was positioned at the same height of As atoms in the 
surface layer. Figures |H1 shows the £DOS at the As atom 
in the surface layer next to the Asc a (a) and As atoms 
distant from the Asca (b). The height of the £DOS just 
below the Fermi energy (between Ep — 0.4 eV and Ep) is 
quite different; As atoms in the surface layer away from 
Asca contain electrons more than an As atom next to 
Asca ■ This means that the electrons on the As atoms 
next to Asca m the surface layer transfer into Asc a ■ 

Atoms in the surface layer align in one-dimensional 
zigzag chains, but each chain does not connect with ad- 
jacent chains. Figure Ufa) shows the iso-surface of the 
highest occupied state which is extended over the whole 
system with the major component along the [110] chain 
of surface As atoms next to the [110] chain containing the 
AsQa ■ The large electron density on this chain originates 
mainly in the dangling bond orbitals of the As atoms. 
The much smaller but some density on the Ascais due 
to the dangling bond p-orbital of the AsQ a atom. The 
electron density on the other As atoms in the [110] sur- 
face chain containing the Asc a does not originate in the 
dangling bond but in the As p-orbitals parallel to the 
dangling p-orbital of the Asca ■ In contrast to the highest 
occupied states, the lowest unoccupied state (Fig. EJb)) 
is localized around the Asc a • The iso-surface consists 
mainly of the dangling bonds of the AsQ a and the As 
atoms in the chain containing the Asca ■ In other words, 
the dangling bond orbitals of the As atoms are occupied 
when AsQa is absent in the zigzag chain, and are unoccu- 
pied when AsQa is present in the chain. This means that 
electrons transferring to AsQ a are those occupying the 
dangling bond orbitals of As atoms next to Asc a in the 
surface in case of the perfect lattice. When Asc a located 
in the surface layer, electron transfer does not occur be- 
cause the energy of dangling bond orbital of As is almost 
the same as that of Asc a and the buckling does not sta- 
bilize the state of the electron transfer. Thus, the surface 
buckling on the site of Asc a in the surface layer does not 
take place. 



V. CONCLUSION 

We presented the electronic structures of the gap state 
associated with a single AsQa in GaAs. The wavefunction 
in the stable configuration mainly consists of the s-orbital 
of AsQa and the p-orbitals of the surrounding As atoms. 
These p-orbitals spread around AsQ a with a radial pat- 
tern of aligned of As p-orbitals. On the other hand, in 
the metastable configuration, the major component of 
the wavefunction is the p-orbital of the Asi and the As 
atoms around it. In comparison with the stable config- 
uration, some p-orbitals of the first neighbor As atoms 
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FIG. 8: The local density of states (ffiOS) of (a) the As atom 
in the surface layer next to AsGa in the surface layer, and (b) 
an As atom in the surface layer further away from AsGa • 
The weight of £DOS just below the Fermi energy (between 
Ef — 0.4 eV and -Ef), of (a) is smaller than that of (b). 



changes their heading directions breaking the radial pat- 
tern of the alignment of the As p-orbitals. As a result, a 
small part of the radial pattern in the stable configuration 
remains. These transitions of the orbital and the bonding 
character can be regarded as the origin of the metasta- 
bility. These essential features of electronic structures 
do not change even when the Asc a is located near the 
surface. However, when the Asc a is located just on the 
surface, the surface unbuckles leaving no localized state 
around the top of the valence band. We can attribute 
the defect image in the STM experiments to these wave- 
function, and the disappearance of the satellite peaks of 
the images to the change in the wavefunction accompa- 
nied by the change to the metastable state. Although 
they does not thoroughly support experimental results, 
the results of the present calculations provide simple and 
clear explanation for the origins of the metastability and 
the characteristic defect images in the STM experiments. 
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FIG. 9: The iso-surface of the electron density on the crystal 
surface of the highest occupied band calculated by the LMTO 
method when an AsGais in the surface layer, (a) Between the 
top of valence band and 0.2 eV below, and (b) between the 
bottom of conduction band and 0.2 eV above it. The AsGa is 
located at the center of the right zigzag chain. 
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